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Preface 


The  enclosed  report  is  a  copy  of  a  paper  submitted  for  publication 
consideration  in  the  Journal  of  Composite  Materials.  The  issuing  of  the 
paper  in  this  report  format  is  for  the  purposes  of  reporting  research 
findings  to  the  agency  which  sponsored  the  work.  The  issuing  of  the  re¬ 
port  is  not  meant  to  duplicate  the  efforts  of  the  Journal  of  Composite 
Material s. 


Abstract 


A  previous  approximate  theory  for  predicting  the  room-temperature 
shapes  of  unsymmetric  laminates  is  examined  in  light  of  the  assumptions 
regarding  the  inplane  strains.  The  previous  theory,  which  was  a  geomet¬ 
rically  nonlinear  extension  of  classical  lamination  theory,  was  felt  to 
be  restrictive  and  this  paper  develops  a  new  theory  in  which  these  re¬ 
strictions  are  relaxed.  It  is  shown  that  despite  the  previous  concern, 
there  is  little  difference  between  the  previous  theory  and  this  theory. 
This  paper  presents  numerical  results  for  the  inplane  residual  strains 
of  unsymmetric  laminates  which  have  cooled  from  curing  into  a  cylindric¬ 
al  room-temperature  shape.  It  is  shown  that  the  residual  strains  are 
compressive  and  practically  independent  of  spatial  location  on  the  lam¬ 
inate.  In  another  facet  of  the  paper,  the  room-temperature  shapes  of 
all  four-layer  unsymmetric  cross-ply  laminates  are  predicted.  There  are 
only  four  unique  stacking  arrangements  for  this  category  of  laminates 
and  it  is  shown  that  their  room-temperature  shapes  are  a  strong  function 
of  their  size  and  their  stacking  arrangement.  Depending  on  these  para¬ 
meters,  the  room-temperature  shape  of  a  four-layer  cross-ply  unsymmetric 
laminate  can  be  a  unique  saddle  shape,  a  unique  cylindrical  shape,  or  a 
cylindrical  shape  that  can  be  snapped  through  to  another  cylindrical 
shape. 


Introduction 


Recently  Hyer  reported  on  a  series  of  experiments  [1]  and  a  theory 
[2]  aimed  at  understanding  the  phenomena  which  govern  the  room-tempera¬ 
ture  shapes  of  unsymmetri cal ly  laminated  composites.  Specifically,  the 
studies  were  concerned  with  [C^/sr^y  and  [O^O^]-]-  laminates  which  are 
fabricated  on  a  flat  plate,  cured  at  an  elevated  temperature,  and  al¬ 
lowed  to  cool  to  room  temperature.  The  studies  were  initiated  because 
it  had  been  reported  by  many  investigators  that  the  room-temperature 
shapes  of  unsymmetric  laminates  did  not  always  conform  to  the  predic¬ 
tions  of  classical  lamination  theory.  Instead  of  being  saddle  shaped, 
as  classical  lamination  theory  predicts,  the  room-temperature  shapes  of 
unsymmetrically  laminated  composites  are  often  cylindrical  in  nature. 

In  addition,  a  second  cylindrical  shape  can  sometimes  be  obtained  from 
the  first  by  a  simple  snap-through  action.  The  experiments  quantified 
some  of  these  effects  while  the  theory  was  developed  to  explain  the 
mechanics  of  the  laminates. 

By  extending  classical  lamination  theory  to  include  geometric  non- 
linearities,  many  of  the  qualitative  and  quantitative  characteri sties  of 
unsymmetric  laminates  were  explained.  The  extended  theory  was  approxim¬ 
ate  in  nature,  using  assumed  displacement  functions  in  conjunction  with 
a  Rayleigh-Ritz  minimization  of  the  total  potential  energy.  For  a  given 
lay-up,  the  theory  predicts  [02/902]t  and  [04/904]^  laminates  to  have 
either  a  single  unique  shape  or  one  of  three  possible  shapes.  The  spe¬ 
cific  shape  depends  on  the  size  of  the  laminate.  When  the  sidelength  of 
a  square  laminate  is  smaller  than  a  certain  critical  length,  a  unique 
saddle  shape  is  predicted.  When  the  sidelength  is  larger  than  this 
critical  value,  either  a  saddle  or  one  of  two  possible  cylindrical 
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shapes  is  predicted.  The  thickness  of  the  laminate  also  has  a  bearing  on 
the  saddle-cylinder  behavior.  A  stability  analysis  shows  that  when 
three  shapes  are  predicted,  the  saddle  shape  corresponds  to  an  unstable 
equilibrium  configuration.  Both  of  the  two  cylindrical  shapes  of  the 
triple-valued  solution  correspond  to  stable  equilibrium  configurations. 
When  the  saddle  shape  is  the  only  shape  predicted,  it  is  found  to  be 
stable.  Correlation  between  the  theory  and  the  limited  amount  of  exper¬ 
imental  data  is  good.  Recently  Hahn  [3]  and  Hahn  and  Hwang  [4]  have 
reported  on  further  experiments  designed  to  check  the  theory.  They 
studied  the  effect  of  sidelength  on  the  shapes  of  [04/90417  laminates 
and  found  reasonable  agreement  between  Hyer's  theory  and  their  experi¬ 
ments. 

While  the  series  of  experiments  Hyer  conducted  involved  a  wide  var¬ 
iety  of  unsymmetric  lay-ups,  his  theory  was  restricted  to  the  study  of 

[02/9021T  and  [04/904]!  laminates.  This  limitation  was  imposed  in  order 
to  keep  the  algebra  at  a  minimum  but  still  retain  the  basic  idea  of 
asymmetry.  In  arriving  at  the  theoretical  predictions,  Hyer  used  simple 
polynomials  to  approximate  the  three  components  of  displacement  in  the 
Rayl eigh-Ritz  schema.  Use  of  these  polynomials  led  to  two  coupled  non¬ 
linear  algebraic  equations.  These  equations  were  solved  numerically  and 
the  solution  of  these  equations  led  directly  to  the  predicted  shapes. 

The  polynomials  used  in  the  displacement  approximati ons  were  of  a  low 
order  so  that  the  number  of  nonlinear  equations  would  be  kept  at  a  mini¬ 
mum.  Because  of  the  good  correlation  between  the  theory  and  the  avail¬ 
able  experimental  data,  the  polynomials  chosen  apparently  modeled  the 
shape  behavior  correctly.  The  question  was  raised,  however,  about  the 
use  of  more  sophisticated  polynomials,  particularly  for  approximating 
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the  inplane  displacements.  The  out-of-plane  displacements  were 
associated  with  the  shapes  and  experimental  observation  led  to  the  form 
of  the  approximating  functions  for  these  displacements.  It  was  found 
from  the  experiments  that  polynomials  were  sufficient.  However,  the  in¬ 
plane  displacements  were  not  as  observable  and  so  there  was  a  lack  of 
physical  evidence  for  choosing  any  particular  functional  form  for  those 
displacements.  Polynomials  were  chosen,  as  opposed  to  other  functional 
forms,  so  that  all  of  the  approximating  functions  would  at  least  be  of 
similar  algebraic  form  and  thus  ease  some  of  the  mathematical  and  alge¬ 
braic  manipulations.  As  is  well  known,  when  using  a  classical  Rayleigh- 
Ritz  approach,  approximating  functions  must  be  chosen  which  are  valid 
over  the  entire  spatial  domain.  The  approximating  functions  must  be  cho¬ 
sen  to  have  enough  freedom  to  accurately  model  complex  behavior  over  the 
domain.  Complex  behavior  may  well  not  exist  but  if  the  possibility  is 
not  allowed  for  in  the  approximate  functional  form,  complex  behavior  may 
be  overlooked  wher  it  does  occur.  For  elastic  systems,  such  as  lamin¬ 
ates,  not  having  enough  freedom  in  the  approximate  displacement  func¬ 
tions  leads  to  a  system  which  acts  much  stiffer  than  the  system  it  is 
modeling.  Thus,  to  examine  the  effects  on  previous  predictions  of 
choosing  more  sophisticated  functions  for  the  approximate  inplane  dis¬ 
placement  fields,  and  to  generate  new  results  for  unsymmetric  laminates, 
the  analytical  work  has  continued  beyond  the  initial  studies.  This  paper 
reports  on  the  additional  work  and  shows  the  effect  of  including  more 
sophisticated  displacement  functions.  Also,  some  new  and  interesting 
results  for  unsymmetric  laminates  are  presented. 

Specifically,  the  paper  examines  the  room-temperature  shapes  of  all 
square  unsymmetric  cross-ply  laminates  which  can  be  fabricated  from  four 
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layers.  There  are  only  four  unique  unsymmetric  laminates  in  this  cate¬ 
gory,  all  other  laminates  being  obtainable  from  these  four  by  simple  ro¬ 
tations.  The  four  unique  laminates  are:  [0/0/0/90]-]-,  [0/0/90/0]r, 
[0/90/0/90]-]-,  and  [0/0/90/90]-]-.  The  last  stacking  arrangement  was  the 
one  examined  in  the  previous  paper.  In  the  notation  used  here,  the 
left -most  entry  in  the  stacking  sequence  nomenclature  corresponds  to  the 
orientation,  relative  to  the  laminate  x  axis,  of  the  fibers  in  the  lam¬ 
ina  at  the  negative-most  z  location.  As  will  be  seen,  [0/90/0/90]-]-  lam¬ 
inates  behave  in  a  manner  similar  to  [0/0/90/90]-]-  laminates.  However, 
[0/0/0/90]t  and  [0/0/90/0]y  laminates  behave  quite  differently.  The  an¬ 
alytical  procedure  used  in  these  latest  studies  parallels  the  approach 
in  the  previous  paper.  For  this  reason  some  of  the  details  of  the  an¬ 
alysis  are  omitted.  This  paper  emphasizes  the  assumptions  used  in  the 
analyses.  The  interested  reader  is  encouraged  to  read  the  previous 
work. 
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Problem  Statement 


Figure  1  illustrates  the  problem  in  question  and  defines  the  co¬ 
ordinate  system  used  in  the  study.  Figure  la  illustrates  the  flat  1am- 
*  iui.e  in  the  uncured  state,  or  in  the  cured  state  at  the  elevated  curing 
temperature.  As  the  laminate  cools  from  the  curing  temperature,  the 
through-the-thickness  lack  of  symmetries  in  the  thermoelastic  properties 
of  the  laminate  causes  it  to  warp  out  of  plane.  The  out-of-plane  dis¬ 
placements  are  accompanied  by  inplane  displacements  but,  as  stated  be¬ 
fore,  these  are  not  as  evident  to  an  observer.  Figure  lb  shows  the 
cooled  laminate  in  a  saddle  shape  and  figs,  lc  and  Id  show  the  laminate 
in  cylindrical  shapes.  As  can  be  seen  in  the  figure,  the  curvatures  of 
the  saddle  are  not  necessarily  equal  nor  are  the  curvatures  of  the  two 
cylinders  necessarily  equal.  This  paper  investigates  the  parameters 
which  influence  whether  the  flat  laminate  of  fig.  la  cures  and  cools  to 
the  shape  given  by  fig.  lb,  lc,  or  Id.  The  solution  procedure  seeks  the 
shape  (or  shapes)  that  minimizes  the  total  potential  energy  of  the  1am- 
i nate. 
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Assimed  Displacement  Functions 


As  in  the  previous  work,  the  approximation  to  the  out-of-plane  dis- 
r'acement,  w,  is  chosen  to  be 


w(x,y)  =  \  (ax2  +  by2). 


(1) 


with  a  and  b  to-be-determined  constants  which  dictate  the  shape  of  the 
laminate.  The  case  a  >  0,  b  <  0  corresponds  to  the  saddle  of  fig.  lb 
while  the  case  a  >  0,  b  =  0  corresponds  to  the  cylinder  of  fig.  lc.  The 
case  a  =  0,  b  <  0  is  another  cylinder  and  is  illustrated  in  fig.  Id.  In 
the  earlier  work,  the  inplane  displacements  of  the  laminate's  midplane 
in  the  x  and  y  directions,  u°(x,y)  and  v°(x,y)  respectively  ,  were  ap¬ 
proximated  using  the  following  logic:  Due  to  the  fact  that  the  out-of¬ 
plane  displacements  of  unsynmetric  laminates  are  many  times  the  laminate 
thickness,  the  von  Karman  approximation  to  Green's  strain  measures 
should  be  used.  Accordingly,  the  strai  n-di  spl  acement  relations  are 
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-  e°x 

z  d2w 

(2) 

5x2 

= 

z  a2w 

(3) 

y 

by2 

=  e° 
bxy 

S2w 

~  Z  dydx  ’ 

(4) 

where  1  is  associated  with  the  x  direction  and  2  is  associated  with  the 

y  direction.  The  inplane  strains,  e°,  e°,  and  e°  ,  are  given  by 

x  y  xy 
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Oil3  I 
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aw  x 

ax  1 


(5) 


3v° 

=  ay  + 


_1  r  SW  <.  2 

2  Vavj 


(5) 


n  l  ,au°  av 0  aw  aw , 

xy  2  ‘  oy  5x  dx  ay 


(7) 


It  '*'3 s  assumed  that  the  vast  majority  of  the  u°  displacement  would  be 
due  to  the  curling  up  of  the  laminate  in  the  x-z  plane.  Referring  to 
fig.  2,  even  though  there  may  not  be  much  strain  in  line  elements  origi¬ 
nally  in  the  x  direction,  these  line  elements  undergo  large  negative 
di sp1 acements  in  the  x  direction.  When  the  laminate  is  flat,  the  dis¬ 
tance  to  point  A,  from  the  origin  0,  is  x.  Vhen  the  laminate  curls  into 
a  cylinder  of  radius  p  with  generators  perpendi cul ar  to  the  x-z  plane, 
the  arc  length  from  0  to  A  is  still  very  close  to  being  x.  However, 
point  A  has  moved  in  the  x-direction  an  amount 


u°  =  psin9-x, 


(81 


where  9  is  the  angular  location,  as  measured  from  the  vertical  z-av* 
of  point  x.  For  cylinders  with  small  angular  openings 


sin  9 


(9) 


By  geometry,  assuming  the  inplane  (arc-wise)  strains  are  negligible,  the 
length  from  0  to  x  is  given  by 


From  en.  1, 


x  =  p9  . 


(10) 
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(11) 


Using  eqs.  9,  10,  and  11  in  eq.  8  leads  to  an  approximate  expres¬ 
sion  for  the  x-direction  displacement  at  point  A  due  to  curling,  namely 


u 


o  = 


(12) 


A  similar  argument  was  made  for  curling  in  the  y-z  plane.  Thus,  the  y- 
direction  displacement  due  to  cylindrical  curling  in  the  y-z  plane  was 
given  by 

v3  b2 

v°(x,y)  =  -  .  (13) 

Continuing  to  construct  u°  and  v°  on  physical  grounds,  it  was  fur¬ 
ther  argued  that  the  inplane  strains  were  similar  to  the  predictions  of 
classical  lamination  theory  in  that  they  were  independent  of  x  and  y. 
Thus  terms  were  added  to  eqs.  12  and  13  which  led  to  constant  inplane 

strains  e°  and  e°.  These  added  terms  were  cx  for  eq.  12  and  dx  for  eq 
x  y 

13,  c  and  d  to-be-determined  constants.  At  this  point  the  inplane 
strains  were  of  the  form 


E 


o 

x 


=  c 


(14) 


e 


0 


y 


=  d 


(15) 


abxy 

• 


(16) 
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As  a  final  step,  it  was  reasoned  that  for  cross-ply  laminates  there 
would  be  no  inplane  shearing  strains.  This  is  the  classical  lamination 
theory  prediction  and  it  was  assumed  to  be  true  for  this  large-deflec¬ 
tion  problem.  Equation  16  contradicted  this  assumption  and  so 
u°  and  v°  were  appended  to  result  in  =  D.  The  final  approximate 
forms  for  u°  and  v°  were 


and 


u°(x,y) 


cx 


a  2  x  3  abxy 2 
6  '  4 


v°(x,y) 


abx  2y 
4 


(17) 


(18) 


where  the  3rd  term  in  each  equation  was  added  to  result  in  zero  inplane 
shear  strain.  These  equations,  along  with  eq.  1,  led  directly  to  the 
following  form  for  the  inplane  strains, 


€ 


0 

x 


aby2 

4 


(19) 


=  d  - 


abx  2 


(20) 


e°  =  0  . 

xy 


(21) 


As  can  be  seen,  e°  could  only  vary  in  the  y  direction  and  e°  could  only 

X  y 

vary  in  the  x  direction.  This  was  felt  to  be  restrictive. 

In  the  latest  work,  to  allow  more  freedom  to  the  spatial  variation 
of  strain,  the  following  functional  forms  have  been  used  to  approximate 
for  the  inplane  elongation  strains, 
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(22) 


= 


du°  1  ,5w 


x  ax 


+  i  ©2  =  al  +  a2x 2  +  a3^2  +  a 4X 2 


o  _  dv°  1 

£y  =  ay  +  2  'ay 


+  2  ^2  =  bi  +  b2*2  +  b3x 2  +  V2*2  •  (23) 


The  a-j  and  are  to-be-determined  constants.  These  functions  are  less 
restrictive  than  the  previously  approximate  strain  fields.  In  the 
present  work  the  elongation  strains  are  allowed  to  vary  in  both  the  x 
and  y  directions.  Only  even  powers  of  x  and  y  have  been  used  in  the 
approximate  strain  functions.  This  is  because  it  has  been  assumed  that 
in  moving  away  from  the  central  portion  of  the  laminate  (the  origin  in 
fig.  la),  the  elongation  strains  increase,  or  decrease,  as  the  edges  of 
the  laminate  are  approached,  independent  of  which  edge. 

Using  w(x,y)  from  eq.  1  in  eqs.  22  and  23  and  integrating  eq.  22 
with  respect  to  x  and  eq.  23  with  respect  to  y  results  in 


ax3  a4x3y2 

u°(x,y)  =  a  jX  +  — —  +  a3xy2  + - 3 —  + 

by 3  b  ^x  2y 3 

v°(x,y)  =  bjy  +  —y  +  b3x2y  + - —  + 


g(y)  - 


h(x) 


bJn3. 

6  * 


(24) 

(25) 


In  these  results  g(y)  is  an  arbitrary  function  of  y  and  h(x)  is  an  arbi¬ 
trary  function  of  x.  For  cross-ply  laminates,  it  is  expected  that 
u°(o,y)  =  0  and  v°(x,0)  =  0.  For  this  to  be  true  it  is  required  that 


g (y )  =  h(x)  =  0  . 

Substituting  eqs.  24  and  25  into  eq.  7  leads  to 


(26) 
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e°y  =  |  { 2( a 3  +  b3)  xy  +  |  (a4x3y  +  tyy3)  +  abxy}  .  (27) 

If,  as  in  the  previous  work,  it  is  assumed  that  the  inplane  shearing 
strains  are  negligible  when  compared  to  the  inplane  elongation  strains, 
then  it  is  necessary  that 

2(a3  +  b3)  +  ab  =  0  ,  (28) 

and 

a4  .  0  »  b4  ,  (29) 

If  it  is  assumed  that 

a 3  -  b^  =  -  4  ,  (30) 


then,  for  the  purposes  of  comparison,  the  present  displacement  field  can 
be  reduced  to  that  of  the  previous  study,  eqs.  17  and  18,  by  setting  a2 
and  b2  bo  zero.  In  addition,  eq.  28  is  also  satisfied.  Using  eq.  30, 
the  assumed  inplane  strain  fields  become 


(33) 


In  the  latest  formulation,  then,  the  inplane  elongation  strains  can  vary 
with  both  the  x  and  y  coordinate.  These  approximate  expressions  for  in¬ 
plane  strains  are  used  in  the  expression  for  total  potential  energy  of 
the  laminate.  The  potential  energy  expression  is  then  subject  to  mini¬ 
mization  with  respect  to  a,  b,  alf  t>i,  a2,  and  b2. 
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Minimization  of  Potential  Energy 


The  total  potential  energy,  W,  for  the  laminate  is  [5] 


W  =  /y0i  “  d Vo  1 ,  (34) 

where  w  is  the  strain  energy  density  and  is  given  by 


u  ~  2  ci  jk1ei  jekl  "  pijeijA^"  * 


(35) 


The  c -j j i  are  the  elastic  constants  of  the  material  and  the  are  con¬ 
stants  related  to  the  elastic  constants  and  the  coefficients  of  thermal 
expansion  of  the  material.  The  temperature  change,  from  curing  to  room 
temperature,  is  aT.  Cooling  corresponds  to  a  negative  AT.  The  ei j  are 
of  course  the  strains.  All  material  properties  are  assumed  to  be  inde¬ 
pendent  of  temperature.  This  is  not  a  restriction  on  the  approach  be¬ 
cause  temperature-dependent  properties  could  be  easily  incorporated  into 
the  theory,  since  aT  *  AT(x,y).  With  the  form  of  eq.  35,  the  datum 
(u)  =  0)  for  the  potential  energy  is  the  elevated  cure  temperature  and 
the  net  work  of  any  external  tractions  acting  on  the  laminate  as  it 
cools  and  deforms  to  its  final  shape  is  assumed  to  be  zero. 

As  is  usually  the  case  with  laminate  analysis,  the  z-direction 
stresses  are  assumed  negligible  and  the  elastic  properties  of  the  indi¬ 
vidual  lamina  are  represented  by  the  familiar  reduced  stiffnesses  re¬ 
ferred  to  the  x-y  system,  0...  The  b..  can  then  be  related  to  the 

ij  ij 

n  's  and  the  lamina  coefficients  of  thermal  expansion  in  the  x  and  y 
i  j 

directions,  a  and  a  respectively.  The  expression  for  the  strain  en- 

x  y 

ergy  density  then  becomes 
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(36) 


co 


=  \  (Qnen2  +  012  ene22  +  2^lle122  +  1  ^22e22 


"  (Oll-x  +  ^12ay^  eil^  ”  (^12°S<  +  ^22ay ^  e22^’ 

Because  these  are  cross-ply  laminates,  there  is  no  a  .  Using  eqs.  1-4 

xy 

and  31-33  in  eq.  36,  the  expression  for  the  total  potential  energy  of 
the  laminate  is  of  the  form 


W  = 


h 

.7 


-h 

z=  Y 


L 

y 

T~ 


L 

x 

.7~ 


x= 


(a,  b,  a br  a2,  b2,  0^ 

<\,  a  ,  AT,  x,  y,  z)  dxdydz  . 


(37) 

With  the  limits  on  the  integral,  it  is  assumed  that  the  laminate  is  of 
thickness  h  and  has  a  length  of  Lx  in  the  x  direction  and  Ly  in  the  y 
di rection. 

The  integral  of  eq.  37,  due  to  the  multiplying  and  squaring  of  the 
strain  expressions,  is  quite  lengthy.  However,  the  spatial  integrations 
on  these  expressions  are  straightforward.  Integration  of  the  variables 
z  and  Q  leads  to  the  familiar  A,  B  and  0  matrices  of  classical  lamina¬ 
tion  theory.  Integration  of  the  variables  z,  Q..,  a  ,  a  ,  and  aT  leads 

ij  x  y 

to  effective  inplane  resultant  forces  and  effective  thermal  moments. 
These  are  the  same  expressions  as  obtained  in  classical  lamination 
theory.  Integration  with  respect  to  variables  x  and  y  results  in  terms 
which  render  the  problem  nonlinear  and  which  dictates  how  the  size  of 
the  laminate  affects  its  room-temperature  shape.  Ultimate  interest  is 
in  the  first  variation  of  W  with  respect  to  the  unknown  variables  a,  b, 
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a  1 ,  b},  a2,  and  b2*  The  theorem  of  minimum  total  potential  energy 
requires  that 


*  •  (f ) 


6a 


*  (f )  a  *  (f 


aa 


1 


6a^  + 


-aw 

■ab. 


+ 


!  aw 


5  o  . 


(38) 


The  actual  taking  of  the  first  variation  results  in  an  equation  of  the 
form, 

6W  =  f^a,  b,  a^  b^  a?,  b?)  5a  +  f2(a,  b,  a^  b  ,  a  ,  b  )  5b 

+  f3(a,  b,  a^  bj,  a2,  b?)  Sa,^  +  (a,  b,  a?,  b?)  fib : 

+  f5(a'  b-  V  Y  a2 ’  V  Y  *  f6  (a'  b-  V  V  V  V  «Y 

(39) 

The  f-j  are  algebraic  relations  between  the  six  unknown  constants  a,  b, 
al,  b i ,  a2,  and  b2,  the  geometric  and  elastic  properties  of  the  lamin¬ 
ate,  and  the  temperature  change  AT.  For  6W  =  0,  the  six  fi  must  all  be 
zero.  The  values  of  a,  b,...etc.  which  make  all  six  f-j  equal  to  zero 
describe  the  equilibrium  shapes  of  the  laminates.  The  functions  f-j  are 
written  out  explicitly  in  the  appendix. 
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Solution  Procedure  and  Stability  Considerations 

Close  examination  of  the  f-j  shows  that  four  of  the  six  of  the  equa¬ 
tions  are  linear  in  a^,  b^,  a^,  and  o2  (these  are  the  first  four  fj  in 
the  appendix).  Using  these  four  equations,  a^,  b^,  a2»  and  b2  are 
solved  for  in  terms  of  a  and  b  and  substituted  into  the  remaining  two 
equations.  This  results  in  two  nonlinear  equations  in  a  and  b.  By 
specifying  the  geometric  and  material  properties  of  a  particular  lamin¬ 
ate,  numerical  values  of  a  and  b  can  be  determined  using  a  Newton-type 
scheme.  The  numerical  values  of  a  and  b  can  then  be  back-substituted  to 
obtain  values  for  a-^  b i ,  a2,  and  b2. 

As  with  the  earlier  study,  there  are  either  three  real  solutions  or 
one  real  solution  to  the  set  of  equations.  Each  solution  corresponds  to 
a  different  room-temperature,  or  equil  ibriun,  shape.  To  complete  the 
solution  process,  each  equilibrium  solution  has  to  be  checked  for  sta¬ 
bility.  This  is  a  standard  procedure  for  any  nonlinear  analysis  which 
involves  multiple  equilibrium  configurations.  The  form  of  the  problem 
here,  eq.  39,  is  of  the  form  discussed  by  Simitses  [6].  For  an  equilib¬ 
rium  configuration  to  be  stable,  the  following  matrix  has  to  be  positive 
def i ni te 


dfi 

afl 

afi 

6f, 

6f, 
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da. 

db. 

da- 

db- 
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da 
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6bj 

da  2 
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6f3 

6f3 
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da 
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dbi 

da^ 

*2 

(40) 
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Using  the  definitions  of  the  fj,  explicit  algebraic  expressions  are  ob¬ 
tained  for  each  term  in  the  matrix.  For  each  equilibrium  solution,  the 
numerical  values  of  a,  b,  a^,  bj,  and  b2  are  substituted  into  each 
term  of  the  matrix  and  positive  definiteness  assessed.  If  the  matrix  is 
positive  definite,  the  equilibrium  configuration  is  stable. 


Numerical  Results 


In  what  follows,  all  numerical  results  are  based  on  laminates  with 
the  following  material  properties: 

El  =  181  GPa  (26.2  x  106  psi ) 

E2  =  10.3  GPa  (1.49  x  10«  psi) 

2  =  0.28 

=  7.17  GPa  (1.04  x  106  Psi) 
a  =  -0.106  x  10-6/°C/  (-0.059  x  10-6/°F) 

=  25.6  x  10"6/°C  (14.2  x  10-6/°F)  . 

The  sidelengths  Lx  and  Ly  are  assumed  to  be  equal  (square  lamin¬ 
ates)  so  that  Lx  =  Ly  =  L.  The  shapes  are  studied  as  a  function  of 
stacking  sequence  and  sidelength.  The  sidelengths  are  in  the  range 
0  <  L  <  800  mm  ( =  32  in). 

Figure  3  shows  the  shape  characteristics  of  square  lamin¬ 

ates  as  a  function  of  sidelength.  Note  the  break  in  the  horizontal 
scale  on  this  figure.  The  upper  portion  of  fig.  3  shows  the  x-direction 
curvature,  a,  as  a  function  of  L  the  lower  portion  shows  the  y-direction 
curvature,  b,  as  a  function  of  L.  Since  these  are  cross-ply  laminates, 
a  and  b  are  the  principal  curvatures.  In  the  x-y  coordinate  system 
there  is  no  twist  curvature.  Also  note  that  the  relation  between  b  and 
L  is  identical  to  the  relation  between  a  and  L  except  for  sign.  The 
solution  characteristics  have  several  branches.  These  branches  are 
identified  by  the  letters  A,  B,  C,  0,  and  E.  It  should  be  noted  that  as 
L  increases,  the  curvature  values  of  branch  BD  approach  zero.  In  addi¬ 
tion,  for  L  >  300  mm  (=  12  in.)  the  curvature  values  for  branch  BE  of 
the  upper  figure  and  branch  BC  of  the  lower  figure  are  actually  zero. 

In  the  figures  these  branches  have  been  separated  from  the  horizontal 
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axes  to  show  their  existence.  The  horizontal  axes  have  been  broken  to 
show  the  details  of  the  curvature  characteristics  for  sidelengths  less 
than  100  nn  (s  4  in.).  As  a  result,  the  distinction  between  branch  BD 
and  the  zero-curvature  branches  becomes  vague  "ear  the  right  end  of  the 
figures.  However,  the  branches  are  distinct  for  all  L.  It  should  be 
pointed  out  that  the  solution  corresponding  to  L  =  9  is  the  classical 
lamination  theory  solution. 

The  shape  character!- sti cs  shown  in  fig.  3  are  virtually  indistin¬ 
guishable  from  the  shape  characteri sties  obtained  using  the  more  limited 
approximate  displacement  fields,  eqs.  17  and  IB.  Thus,  the  more  sophis¬ 
ticated  inplane  functional  forms  do  not  affect  the  predicted  out-of¬ 
plane  displacements  to  any  significant  degree. 

The  stability  analysis  shows  branch  AB  represents  a  stable  equilib¬ 
rium  solution.  Since  b  =  -a  on  this  branch,  the  branch  represents  lam¬ 
inates  which  exhibit  saddle  shapes  at  room  temperature.  Branch  BD,  a 
continuation  of  the  saddle  branch,  is  shown  to  be  unstable.  Branch  BC 
is  stable  and  corresponds  to  laminates  which  have  a  large  curvature  in 
the  x  direction  and  a  small  curvature  in  the  y  direction.  As  the  side- 
length  L  increases,  the  y-direction  curvature  asymptoti cal ly  approaches 
zero.  This  branch  represents  circular  cylinders  with  generators  in  the 
y  direction.  Branch  BE  is  also  stable  and,  in  the  limit  of  large  L, 
corresponds  to  laminates  with  large  curvature  in  the  y  direction  and  no 
curvature  in  the  x  direction,  i.e.,  cylinders  with  generators  in  the  x 
direction.  Thus,  a  square  [02/902]!  laminate  with  a  sidelength  less 
than  35  mm  (~  1.4  in.)  will  cure  to  the  saddle  shape.  Generally,  for 
non-zero  sidelength  the  saddle  will  be  shallower  than  the  saddle  predic¬ 
ted  by  classical  lamination  theory.  For  laminates  greater  than  75 


18 


mm  (=  3  in.)  on  a  side,  the  room-temperature  shape  will  be  one  of  two 
cylindrical  shapes.  To  get  from  one  shape  to  the  other  requires  a  snap- 
through  action.  Laminates  with  sidelengths  in  the  range  35-75  mm  will 
have  two  possible  shapes  and  a  snap-through  action  is  also  required  to 
get  from  one  shape  to  another.  Neither  shape  will  have  associated  with 
it  a  direction  of  zero  curvature.  For  either  shape,  one  curvature  will 
be  much  smaller  than  the  other  and  as  the  laminate  sidelength  increases 
beyond  75  mm,  the  smaller  curvature  will  asymptotically  approach  zero. 

Although  the  effect  of  the  more  sophisticated  inplane  displacement 
fields  On  the  predicted  out-of-plane  deformations  is  negligible,  their 
effects  on  inplane  strains  could  be  important.  Surpri si ngly,  the  effect 
is  also  negligible.  Figure  4  shows  the  predicted  inplane  residual 
strains,  and  as  a  function  of  spatial  location  for  two 
laminates.  One  laminate  has  a  sidelength  of  150  mm  ( «  5  in.)  while  the 
other  laminate  has  a  sidelength  of  300  rm  ( «  32  in.).  Based  on  the  re¬ 
sults  of  fig.  3,  both  laminates  are  cylindrical.  For  the  cases  shown, 
the  cylinders  have  generators  in  the  y  direction  (fig.  lc  and  branch  BC 
of  fig.  3).  The  strains  were  computed  using  eqs.  31  and  32,  functional 
forms  which  allow  spatial  variations  of  strain  in  both  the  x  and  y  dir¬ 
ection.  As  can  be  seen,  for  the  larger  cylinder  both  inplane  elongation 
strains  are  spatially  invariant.  For  the  smaller  laminate,  varies 


slightly  in  the  y  direction  and  is  independent  of  x.  On  the  other  hand, 
for  the  smaller  laminate,  varies  slightly  with  x  but  is  independent 
of  y.  Fxamination  of  eqs.  10  and  20,  tho  strains  assumed  in  the 
previous  theory,  shows  that  neither  the  x  variation  of  e?  nor  the  y  var¬ 
iation  of  was  allowed.  This  apparently  was  not  as  restrictive  as 
first  felt  since,  from  fig.  A,  there  is  no  tendency  of  e°  to  vary  with  x 
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or  to  vary  with  y.  In  fact,  the  strains  predicted  using  the  previous 

theory  are  indistinguishable  from  the  strains  predicted  by  the  present 

theory.  As  a  matter  of  interest,  classical  lamination  theory  predicts 

e°  and  e°  to  be  -905  x  10"^  for  these  cases.  The  elongation  strains  as- 
x  y 

sociated  with  the  direction  of  nonzero  curvature  are  very  close  to  the 
predictions  of  classical  lamination  theory.  The  strains  associated  with 
the  direction  of  zero  curvature  are  lower.  As  will  be  pointed  out 
later,  the  curvature  of  the  cylinder  is  also  roughly  the  value  of  one 
curvature  of  the  saddle  predicted  by  the  classical  theory.  This  can  be 
seen  by  noting  in  the  upper  portion  of  fig.  3  that  branch  AC  asymptotes 
close  to  the  value  given  by  point  A  (L  =  0),  the  classical  theory  curva¬ 
ture  prediction.  The  roles  of  c°  and  e°  are  reversed  for  the  other  cyl  - 

x  y 

indrical  configurations  at  these  lengths,  i.e.  a  =  0  b  <  0. 

From  the  results  of  fig.  4,  and  from  other  strain  studies  not  shown 
here,  it  appears  that  the  larger  the  cylindrical  laminate,  the  more  spa¬ 
tially  uniform  the  inplane  strains  become,  it  is  important  to  note  that 
the  inplane  strains  are  everywhere  compressive.  The  laminate  not  only 
has  a  useful  shape,  a  cylinder,  but  it  is  prestrained  in  compression. 
This  compressive  prestrain  would  be  very  useful  if  the  cylinder  were 
used,  for  example,  as  an  internally  pressurized  vessel.  Owing  to  the 
curvature,  the  strains  away  from  the  midsurface  vary  with  through-the- 
thickness  location.  However,  foe  the  cases  in  fig.  4,  all  points  in  the 
laminates  are  in  a  state  of  compressive  residual  strain.  The  largest 
compressive  strains  occur  in  the  matrix  direction  at  the  inner  surface 
of  the  cylinder.  The  smallest  compressive  strains  occur  in  the  fiber 


direction  at  the  outer  surface.  For  thicker  laminates,  a  [04/904]]-  for 
example,  tne  thickness  and  curvature  effects  combine  to  produce  tensile 


strains  in  the  fibers  at  the  outer  surface. 


Figures  5-7  show  the  shape  characteristics  of  the  other  unsymmetric 
cross-ply  laminates  studied.  All  laminates  are  also  square  with  the 
curvature  characteristics  in  the  x  and  y  direction  being  illustrated  as 
a  function  of  sidelength.  There  is  also  a  break  in  the  horizontal  scale 
of  fig.  5.  The  horizontal  scales  of  figs.  6  and  7  are  continuous  but 
the  vertical  scales  of  bran  ;h  DBC  in  each  figure  have  been  exaggerated. 
For  all  laminates  a  and  b  are  the  principal  curvatures.  Stability  con¬ 
siderations  shows  branch  BD  to  be  unstable  for  all  laminates. 

The  [0/90/0/90]t  shape  cha^acteri sties  are  similar  to  the  shape 
characteri sti cs  of  the  [0/0/90/90] -p.  However,  the  shape  characteristics 
of  the  [0/0/0/90] j  and  the  [0/0/90/0] 7  are  quite  different.  The  most 
striking  difference  is  the  lack  of  a  coalescence  of  the  solution  bran¬ 
ches.  This  coalescence  of  solutions,  which  is  denoted  as  point  8  in 
figs.  3  and  5,  is  often  referred  to  as  a  bifurcation  point.  The  dis¬ 
joint  behavior  of  the  solution  branches  for  the  [0/0/0/90] T  anq 
[0/0/90/  0]t  1S  °ften  referred  to  a  limit  point  behavior.  Generally,  the 
existence  of  a  bifurcation  point,  as  opposed  to  the  existence  of  a  limit 
point,  is  associated  with  the  analysis  of  a  perfect  or  ideal  case.  In 
the  problems  here,  the  perfect  cases  are  associated  with  'perfect  unsym¬ 
metry'  or  'antisymmetry.'  The  [0/0/90/90]t  and  [0/90/0/90] j  are  anti¬ 
symmetric.  In  contrast  a  [l/0/90/90]y,  for  example,  would  exhibit  limit 
point  behavior  rather  bifurcation.  The  perfect  asymmetry,  or  the  lack 
of  it,  is  reflected  in  ie  A,  B,  0  matrices.  Whereas  as  [0/0/90/90] T 
laminates  cease  to  have  a  unique  stable  saddle  shape  when  the  sidelength 
exceeds  35  mm,  [0/90/0/90]-]-  laminates  retain  a  unique  saddle  shape  up  to 
sidelengths  of  about  90  mm  (=  3.5  in.).  This  is  a  reflection  of  the 
fact  that  the  [0/0/90/90]j  configuration  has  a  greater,  or  more  serious. 
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asymmetry  than  does  the  [0/90/0/90]y  lay-up.  The  level  of  asymmetry  is 
reflected  in  the  8  matrix.  For  sidelengths  greater  than  200  mm 
«  8  in.)  branches  BC  and  BE  of  fig.  b  represent,  respectively,  cylinders 
with  generators  in  the  y  and  x  directions.  To  get  from  one  cylindrical 
shape  to  another  requires  a  snap-through  action.  Like  the  [0/0/90/90]j 
cases,  the  relation  between  b  and  L  is  identical  to  the  relation  between 
a  and  L  except  for  sign. 

The  [0/0/0/90]j  and  the  [0/0/90/0]y  are  quite  interesting  and  de¬ 
serve  some  discussion.  Of  course,  at  zero  sidelength  the  classical  lam¬ 
ination  theory  saddle  shape  is  predicted.  Unlike  the  previous  two 
cases,  the  curvatures  of  the  saddle  are  not  equal.  As  sidelength  in¬ 
creases,  the  smaller  of  the  two  curvatures,  a,  is  quickly  suppressed  to 
zero.  The  larger  curvature,  b,  essentially  remains  unaffected  by  the 
sidelength  dimension.  For  sidelengths  between  100  and  180  mm,  ( =■  4  and 
»  7  in.),  [0/0/0/90]]-  laminates  cure  to  a  cylindrical  shape.  What  is  is 
important  is  that  this  cylindrical  shape  is  the  only  shape  that  can 
exist  at  those  sidelengths.  The  cylindrical  shape  is  unique  and  it  is 
stable.  For  the  'perfect'  cases,  i.e.  the  [0/0/90/90]y  and 
[0/90/0/90]t>  at  sidelengths  for  which  the  cylindrical  behavior  appears, 
one  cylinder  can  be  snapped  into  another  cylinder.  For  these  cases  the 
cylindrical  shape  appears  only  in  the  presence  of  another  possible  cyl - 
indrical  shape  and  is  thus  a  bistable  configuration.  From  the  results 
of  fig.  6  it  appears  there  is  the  potential  for  generating  unique  cylin¬ 
drical  shapes  with  unsymmetric  laminates.  Similar  comments  can  be  made 
concerning  the  uniqueness  of  shapes  of  [0/0/90/0]y  laminates  with  side- 
lengths  between  120  and  525  mm  ( =>  5  and  =>  21  in.). 
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For  sidelengths  greater  than  the  length  associated  with  point  B, 
both  the  [0/0/0/90]]-  and  the  [0/0/90/0]]-  laminates  exhibit  two  possible 
cylindrical  shapes.  The  two  cylindrical  shapes  possible  with  each  lam¬ 
inate  have  curvatures  of  opposite  sign  and  unequal  magnitude.  These 
cylindrical  shapes  are  bistable.  It  is  not  clear,  however,  that  the 
amount  of  force  required  to  make  the  snap  from  BC  to  AE  is  the  same  as 
the  force  required  to  make  the  reverse  snap.  The  equality  of  snap- 
through  force  is  most  likely  the  case  for  the  [0/0/90/90]y  and 
[0/90/ 0/903 y  laminates. 

There  are  some  interesting  similarities  among  the  four  laminates 
considered.  In  all  cases  all  branches  are  stable  except  branch  BD.  In 
each  case  BD  represents  all  or  part  of  a  saddle-shaped  configuration. 

It  appears,  then,  the  shape  qualitatively  predicted  by  classical  lamina¬ 
tion  theory  is  often  inherently  an  unstable  configuration.  Also  for  all 
cases,  as  laminate  sidelength  increases  slightly  from  zero,  say  up  to  3D 
mm  (1.2  in.),  the  curvatures  seem  to  decrease.  For  [0/0/0/90] j  and 
[0/0/90/0]y  laminates,  the  smaller  of  the  two  curvatures,  the  x-direc- 
tion  curvature,  decreases  rapidly  toward  zero.  1  me  major  curvature,  the 
y-direction  curvature,  decreases  slightly.  In  the  earlier  experimental 
studies  [1]  this  suppression  of  the  smaller  curvature  was  observed  to 
always  be  the  case.  For  [0/0/90/901]  and  [0/90/0/90]y  laminates,  as 
sidelength  increases  slightly  from  zero,  both  the  x  and  y  curvatures  de¬ 
crease  considerably  but  they  do  not  go  to  zero.  In  all  cases,  as  side- 
lengths  increase  beyond  30  mm,  if  one  of  the  curvatures  is  not  zero,  it 
eventually  becomes  zero.  The  other  curvature  eventually  returns  to  a 
value  close  to  the  value  predicted  by  classical  lamination  theory.  Spe¬ 
cifically,  for  both  the  [0/0/90/90]]-  and  [0/90/0/90]]-  laminates,  the  y- 
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direction  curvature  decreases,  in  magnitude,  and  then  increases  along 
path  ABE,  as  in  the  lower  portion  of  fig.  3.  The  x-direction  curvature 
decreases  and  is  ultimately  suppressed  to  zero  along  path  ABE,  as  in  the 
upper  portion  of  fig.  3.  For  both  the  [0/0/0/90]j  and  the  [0/0/90/0]y 
laminates,  the  y-direction  or  major  curvature  decreases  slightly  in  mag¬ 
nitude  and  then  increases  along  path  AE,  as  in  the  lower  portion  of  fig. 
6.  The  x-direction  or  minor  curvature  is  suppressed  to  zero  along  AE, 
as  in  the  upper  portion  of  fig  6.  For  the  perfectly  asymmetric  cases, 
path  ABC  exhibits  similar  decreasi ng-suppressi on  and  decreasi ng-i ncreas- 
ing  curvature  behavior  for  the  two  curvatures  involved. 

Though  not  shown,  the  residual  strain  levels  for  all  four  laminates 
are  similar.  There  is  nothing  unusual  about  the  spatial  variation  of 
strain  for  any  of  the  laminates  considered.  Although  the  magnitudes  of 
the  strains  are  different  among  the  four  laminates  studied,  fig.  3  is 
representative  of  the  strain  behavior  for  all  cylindrical  configurations 
of  all  laminates. 


24 


Acknowl edgements 


The  work  reported  on  here  was  conducted  while  the  author  was  being 
supported  by  NASA  Cooperative  Agreement  NCCI-15  and  AFOSR  Grant  Rl-0195. 
The  financial  support  of  the  two  funding  agencies  is  gratefully  acknow- 
1  edged. 


25 


References 


1.  Hyer,  M.  W. ,  "Some  Observations  on  the  Room-Temperature  Shapes  of 
Unsymmetrical 1y  Laminated  Composites,"  J.  Comp.  Materials,  Vol .  15, 
pp.  175-194,  1981. 

2.  Hyer,  M.  W. ,  "Calculation  of  the  Room-Temperature  Shapes  of  Unsym- 
metrically  Laminated  Composites,"  J.  Comp.  Materials,  Vol.  15,  pp. 
296-311,  1981. 

3.  Hahn,  H.  T. ,  "Warping  of  Unsymmetric  Cross-Ply  Graphite/Epoxy  Lamin 
ates,"  Composites  Technology  Review,  Vol.  3,  No.  3,  pp.  114-117, 
1981. 

4.  Hahn,  H.  T.  and  Hwang,  D.  G. ,  "Residual  Stresses  and  Their  Effects 
in  Composite  Laminates,"  Proceeding  of  Symposium  on  Engineering  Sci 
ence  and  Mechanics,  Tainan,  Taiwan,  Dec.  1981. 

5.  Fung,  Y.  C.  Foundations  of  Solid  Mechanics,  Prentice-Hal 1 ,  Inc., 
Englewood  Cliffs,  Nj,  1965,  p.  354. 

6.  Simitses,  G.  J.  An  Introduction  to  the  Elastic  Stability  of  Struc¬ 
tures,  Prentice-HaTl 7  IhcTTngl ewood  Cliffs,  NJ,'  1976,  pp.  8-147 


APPENDIX 


GOVERNING  ALGEBRAIC  EQUATIONS 

The  following  six  equations  govern  the  equilibrium  shapes  of  the  lamin¬ 
ate. 


fl  =  Vl  +  V2  +  A12»l  +  4C6b2 


A-l 


-  Nx  -  C^ab  -  B^a  -  C^ab 


f2  '  Vl  +  4C4a 2  +  A22bl  *  S(lb2 


A-2 


-  N  -  Cfiab  -  C7ab  -  82,5b 


VSV'mWi  *  TWC5b2 

'  \T  -4"  '  Glat>  "  G2a  *  ®6ab 


A-3 


f 4  =  4C6al  +  144  C5a2  +  G10bl  +  G15b2 

T  L  2 

-  N  — -  G3ab  -  Gyab  -  Ggb 


A-a 


f5  =  -  C^b  -  Gla2b  +  C2ab2  ♦  C3ab  -  Buai 


"  G232  +  C3ab  +  °Ua  '  CRblb  '  Sb2b  +  C5ab2 


A-5 


+  D12b  -  C4aib  -  Gba2b  +  C53b2  *  G7blb  -  G7b2b 


L  2 


L  2 


*  C„ab2  *  C0t>2  *  (N/  -*glb  *  fb„T  *-„)b  *  II 


'y  4R 
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The  constants  A 


11,  A12,  A22,  B 1 1 ,  B22 »  Dll,  r22>  and  D12  have  their  usu¬ 
al  meaning.  The  NXT,  NyT,  MXT,  and  MyT  are  effective  inplane  thermal 
forces  and  moments,  respectively,  and  are  given  by: 


T  h  /  2  h  /  2 

Nx  =  aT  -h/2  tillV  ^12ay  ’  Ny  =  AT  ■^-h/2^12ax  +  ^22ay^dz 

Mx  =  AT  ^h/2^11ax  =  ^I2ay^zdz:  My  =  aT  ^-h/2^12°x  +  ^22ay^zdz 


29 


Fig.  1.  Laminate  shapes:  (a)  at  the  elevated 
curing  temperature,  and  at  room  temperature, 
(b)  a  saddle  shape,  (c)  a  cylindrical  shape, 
and  (d)  another  cylindrical  shape. 


Fig.  2.  Kinematics  of  a  flat  laminate  deforming 


into  a  cyl inder. 


cure  temperature 


Fig.  3.  Room-temperature  shapes  of  square 
[0/0/90/90]t  laminates. 


Fig.  4.  Spatial  dependence  of  residual  inplane  strains  e°  and  e° 

x  y 

for  two  CG2/9O2DT  cylindrical  laminates. 


Fig.  5.  Room-temperature  shapes  of  square 
[0/90/0/90] -p  laminates. 


Fig.  6.  Room-temperature  shapes  of  square 
[0/0/0/90] j  laminates. 


m 


Fig.  7.  Room  temperature  shapes  of  square 
[0/0/90/0]t  laminates. 


